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Atomic clocks have recently reached a fractional timing precision of < 10“^®. We point out that an 
array of atomic clocks, distributed along the Earth’s orbit around the Sun, will have the sensitivity 
needed to detect the time dilation effect of mHz gravitational waves (GWs), such as those emitted 
by supermassive black hole binaries at cosmological distances. Simultaneous measurement of clock- 
rates at different phases of a passing GW provides an attractive alternative to the interferometric 
detection of temporal variations in distance between test masses separated by less than a GW 
wavelength, currently envisioned for the eLISA mission. 
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Introduction. Over the past year, the precision of op¬ 
tical lattice clocks has advanced dramatically, to a frac¬ 
tional timing precision of (At/t) ^ 10“^®, with prospects 
for a future improvement by two additional orders of 
magnitude through the use of other atoms [in. Here, 
we point out that the new regime of timing precision 
made accessible by atomic clocks overlaps with the ex¬ 
pected amplitudes of time dilation and compression due 
to the passage of gravitational waves (GWs). The stan¬ 
dard time-dilation effect for a clock at some distance from 
a black hole, would be modulated by the periodic change 
in this distance due to the orbital motion in a binary 
black hole system. Quantitatively, compact binaries of 
supermassive black holes ii at cosmological distances 
produce a periodic modulation of the time-time compo¬ 
nent of the metric (in a Newtonian gauge) at a level of 
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where Dl{z) is the luminosity distance for a source at 
redshift z, M^ = (l + /{Mi + Ma)!/® with 

Ml and M 2 being the masses of the binary members 
which are assumed here to be on a circular orbit, and 
/ = {w/2'k) is the observed (redshifted) GW frequency 
(obtaining its maximum value based on the orbital fre¬ 
quency at the binary’s innermost stable circular orbit).^ 
Such black hole binaries are a natural consequence of 
galaxy mergers Q- Around the Milky Way’s super¬ 
massive black hole, GWs with a similar amplitude and 
frequency can be emitted by orbiting close-in stars Q. 
Such sources constitute th^rimary targets for the future 
eLISA space observatory 0, which aims to detect GWs 


^ In this paper, we adopt for pedagogical reasons a Newtonian 
gauge which is commonly used to describe the time-dilation ef¬ 
fect due to stationary gravity, as measured in the Pound-Rebka 
experiment Q. In this gauge, an oscillating perturbation in the 
time-time component of the metric, hoo, would trigger periodic 
variation in the Pound-Rebka time dilation and a mismatch be¬ 
tween the ticking rate of clocks separated apart. 


in the frequency range of 0.1-100 mHz and is planned for 
launch in two decades. eLISA is designed as a laser inter¬ 
ferometer that will record the phase shift introduced by 
a passing GW as the wave induces a change in the space- 
time curvature, and hence a change in light travel time 
between its three test masses, which are separated from 
each other by 10® km (less than the GW wavelength). 
Here, we consider the alternative approach of detecting 
the differential time dilation experienced by clocks lo¬ 
cated at different phases of a passing GW. For this pur¬ 
pose, we propose using an array of atomic clocks. 

Atomic clocks in space. We envision a set of small 
orbiting units (at a minimum two of them) equipped 
with atomic clocks, and a primary spacecraft between 
them, as illustrated in Fig. 1. For simplicity, the clocks 
are assumed to be distributed along the circular orbit 
of the Earth around the Sun at an orbital radius of 
8.3 light-minutes (1AU= 1.5 x 10®km), since this con¬ 
figuration minimizes the kinetic energy requirements for 
launch. The typical distances between units will thus be 
~ 10® km. Additional units will naturally improve the 
sensitivity and directional angular resolution. 

A passing GW characterized by a metric perturbation 
will induce periodic variations (over a pe¬ 
riod 2tt/uj = 2'K/kc) in the ticking rate of each clock 
dg. The amplitude of the timing variation would be 
ihoo (since time is dilated by (1 -I- /iqo)”^/^). For each 
pair of clocks, the relative phase of the periodic varia¬ 
tion = kAr will be dictated by the projection of 
the difference in clock position vectors Ar on the GW 
wavevector k (which defines the GW propagation di¬ 
rection). The largest timing phase difference would be 
realized between clocks separated by half a wavelength, 
iA = lAU(//lmHz)“^, corresponding to a timing phase 
difference At/) = tt. Measurements of the N{N — 1) phase 
differences 2 n) ^ ^ clocks would 

allow to localize the GW source position on the sky as 
—k. 

To communicate the periodic time dilation or compres¬ 
sion signal associated with a GW, the clock in each unit 
should drive an optical-frequency (zziaser 10^® Hz) laser, 
pointed at the primary spacecraft. The required signal- 
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FIG. 1: Schematic illustration of an array of atomic clocks 
in a circular orbit around the Sun, as an observatory for 
low-frequency GWs. Clocks separated by half a wavelength, 
= tt/A: = lAU(//mHz)~^, of a passing GW, would show 
the largest difference in fractional timing, hoo, varying pe¬ 
riodically over the GW period, f~^ = 2tv jkc. We adopt a 
Newtonian gauge for pedagogical reasons, but the measured 
physical effect should be independent of gauge. 

to-noise ratio is not dictated by the precision for distance 
measurements through laser phase shifts, but rather by 
what is needed to maintain an optical phase lock be¬ 
tween the different units such that the ticking rate of their 
clocks can be compared at a high precision. Techniques 
for remote optical clock comparison had been developed 
over the past decade [ll|, [l^ and are conceptually dif¬ 
ferent from the interferometric technique currently envi¬ 
sioned for eLISA. In the absence of phase noise in space, 
the distance noise level of ^ 5 x 10“® cm/-\/Hz expected 
from the shot noise of a 1 Watt laser detected by tele¬ 
scopes of 30 cm diameter across a 1-2AU path length, 
would be more than sufficient to phase lock a single op¬ 
tical cycle. The needed laser is available and already 
in use. The primary spacecraft will interfere the laser 
signals from the two daughter clocks to produce a beat 
signal with a frequency i^heat ^oo^^iaser, be. with a pe¬ 
riod of < 10^ s, comparable to the period of the GW, 
and hence detectable. 

Discussion The proposed method is fundamentally 
different from the current interferometric methods which 
underline the design of Advanced-LIGO [a and eLISA. 
In GW interferometry, the GW wavelength is consid¬ 
erably longer than the interferometer arm length, and 
at any given moment one is essentially measuring the 
roughly-uniform space curvature in the region between 
the test masses by means of the light-travel time. Half a 
GW period later, the change in light-travel time reveals 
the passage of the GW. In our clock-timing method, the 
separation between clocks is of order half a GW wave¬ 
length and one measures the simultaneous difference in 
clock rates between the low-curvature and high-curvature 
phases of the GW (which again flips after half a GW pe¬ 


riod). The change in the light travel time between the 
clocks due to the GW averages out over a GW wave¬ 
length, and thus amounts to a higher order correction to 
the beat frequency, which can be neglected. 

Our proposal is to measure timing variations rather 
than distance variations. Our detection scheme is not 
concerned with laser phase variations that are induced by 
distance variations between free-floating units, but rather 
with the change in the rate at which clocks are ticking 
relative to each other. The lasers are locked to the clocks 
so as to keep their timing stability at a high level. 

The precision achieved by optical lattice clocks scales 
inversely with the square root of the integration time 
0 and limits the clock-timing method to low-frequency 
GWs with / < 1 mHz, allowing integration for ~ ^ 

lO^s. This, in turn, requires clocks at ~ 1 AU separa¬ 
tions, as considered here. Secular trends (due to slow 
orbital variations in the gravitational effect of the Sun or 
planets, the energy loss of the Sun, or the phase shift due 
to the solar wind) over the short period of the measure- 
ment (< 1 hour) can be separated from the periodic GW 
signal in the frequency range of interest here. Further¬ 
more, for GW signals lasting weeks or months, the orbital 
motion will scan the ecliptic, improving the source local¬ 
ization accuracy. Solar oscillations at mHz frequency (p- 
mode and especially g-mode) could also be detected by 
the proposed array. 

As indicated by Eq. o, the timing precision of atomic 
clocks is sufficient to detect a GW from a supermassive 
black hole binary at cosmological distances within one 
GW period. The signal-to-noise ratio (SNR) of GW de¬ 
tection will improve in proportion to the square-root of 
the number of wave cycles being observed. Therefore, to 
achieve SNR= I, the fractional timing precision per cycle 
could be worse than the GW amplitude by the square- 
root of the number of wave cycles observed. 

So far we assumed a monochromatic GW signal and 
ignored the slow drift in the GW frequency (so-called 
chirp) due to the inspiral of the supermassive black hole 
binary. The lifetime of tight binaries is in fact limited 
by the rate of GW energy loss [3. For a circular bi¬ 
nary orbit, the fractional frequency increase during an 
infinitesimal observing time Atobs is given by 



The temporal rise in frequency / (chirp) of the GW sig¬ 
nal can also be detected through clock timing to provide 
a second constraint on the values of Adz and Dl in addi¬ 
tion to Eq. ©■ Since there are > 10^^ galactic mergers 
within a Hubble time in the observable volume of the 
universe, the duty cycle of detectable signals (event rate 
times lifetime) could be high 

The use of a large number of uncorrelated clocks, N, 
in every clock unit, can improve the timing precision by a 
factor of This may become another advantage of 

the proposed method, as atomic clocks become progres¬ 
sively small, low-weight and inexpensive. If the clocks are 
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quantum entangled, the timing precision could improve 
as 1 /N and inversely with integration time m. 
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